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Abstract 

The complexation of xylazine with P-cyclodextrin was studied in aqueous medium. NMR titrations confirmed the formation of a 
1:1 inclusion complex. A ROESY spectrum was recorded with long mixing time which contained TOCSY artifacts. It only 
confirmed the presence of xylazine aromatic ring in the P-cyclodextrin cavity. No information regarding the mode of penetration, 
from the wide or narrow side, could be obtained. We calculated the peak intensity ratio from the inter-proton distances for the most 
stable conformations obtained by molecular dynamics studies in vacuum. The results show that highly accurate structural informa- 
tion can be deduced efficiently by the combined use of quantitative ROESY and molecular dynamics analysis. On the other hand, a 
ROESY spectrum with no spin diffusion can only compliment an averaged ensemble conformation obtained by molecular dynamics 
which is generally considered ambiguous. 



Introduction 

The emergence and establishment of supramolecular chemistry 
as an important domain of science has fueled the development 
of complex chemical systems from components, interacting by 
non-covalent intermolecular forces. This field transcends the 
traditional barriers separating many disciplines of science and is 
the basis for most of the vital biological processes [1]. The basis 
of supramolecular chemistry is molecular recognition where 
host and guest species interact with each other and exist as a 
single system. These host-guest systems symbolize simplest 



examples of supramolecular systems in which a guest is encap- 
sulated into the internal cavity of a larger host molecule. The 
most widely used hosts are cyclodextrins (CDs) which are crys- 
talline, homogeneous and non-hygroscopic substances 
composed of a- 1^4 linked glucose units. The outside surface 
of CDs is hydrophilic while the interior of the cavity is 
hydrophobic [2,3]. The encapsulation of a guest into the CD 
cavity has a profound effect on the chemical, physical and bio- 
logical properties of the guest. 
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The structure establishment of CD inclusion complexes in solu- 
tion state is a challenging task considering the fact that the 
inclusion phenomenon is a dynamic process, the included guest 
being in fast exchange between the free and bound state, and 
hence deriving instantaneous information with regard to struc- 
tural changes is not easy. 

NMR spectroscopy has evolved as a method of choice for 
studying inclusion complexes in solution, although only time- 
averaged structural information can be extracted in NMR time 
scale. Nuclear Overhauser Enhancement (NOE) experiments, 
which depend on intemuclear distances, are used to identify the 
part of guest that is involved in complexation and to determine 
the position of the guest inside the CD cavity. However, the 
inability of NMR spectroscopy lies in the fact that it is not able 
to describe the mode of guest penetration, i.e. through the wide 
or the narrow rim side. The use of computational approaches 
along with NMR studies has proved to hold promise for systems 
which have dynamic nature [4]. The combination of these two 
techniques helps in understanding the complexation process and 
in gaining a deeper insight into the geometry of the system with 
a high degree of accuracy [5] and several cyclodextrin 
complexes have been investigated using this approach [6,7]. It 
has been observed that by the use of molecular dynamic tech- 
niques with no prior knowledge about the existence or geom- 
etry of CD inclusion complexes [8], accurate and reliable 
predictions can be made regarding the possible formation and 
other aspects of inclusion process by following a general simu- 
lation protocol [9]. 



Results and Discussion 
NMR studies 

Demarco and Thakkar first observed the chemical shift changes 
of the H-3' and H-5' protons, positioned inside the cavity of 
P-CD, in the presence of aromatic compounds which is inferred 
due to anisotropic effect of the aromatic ring entering the cavity 
[10]. Numerous studies ever since have shown that cavity 
protons and sometimes H-6' also, located on the cavity rim at 
the narrow end, undergo appreciable shift changes upon inclu- 
sion of the guest while H-1 ', H-2' and H-4', located outside the 
cavity, are relatively unaffected [11-14]. 

The ^H NMR spectral data of P-CD were in good agreement 
with the reported [15,16] and each proton resonance of 
xylazine, especially of the aromatic region, was assigned with 
the help of COSY data. The aromatic protons of xylazine were 
observed as a distorted doublet at 7.13 (2H, J=9 Hz, H-1) and 
a doublet of doublet at 7.20 (IH, = J2 = 9 Hz, H-2). 
Remaining protons were found resonating at 2.10 (H-4, -CHS) 
and 3.25 (H-3, 5). 

^H NMR spectra of mixtures of xylazine and P-CD displayed 
upfield shift changes (A5), compared to pure P-CD, in the H-3 ' 
and H-5' proton resonances which were affected by concentra- 
tion of xylazine. The magnitude of A5h-3' was slightly larger 
than A5h-5'- Signals for H-2' and H-4' also slightly moved 
highfield but were unaffected by the change in concentration of 
xylazine (Figure 2). On the other hand, changes in shape and 
size of chemical shifts were observed for aromatic protons of 



We report here the study of the complexation of P-CD with 
xylazine (XZ) (Figure 1) in aqueous medium by NMR spec- 
troscopy and molecular dynamics. Xylazine has two rings, one 
substituted benzene ring which can only partially penetrate the 
CD cavity due to the presence of two methyl groups and a hete- 
rocyclic ring which should prefer aqueous medium outside the 
cavity. The study suggests that highly reliable structural infor- 
mation can be deduced by a combination of a poorly resolved 
^H-^H ROESY spectrum, which otherwise gives only the infor- 
mation about the part of the guest included, and careful inter- 
pretation of molecular dynamics results in vacuum though it 
does not fully simulate the experimental conditions. 
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Figure 1: Structure of xylazine. 



4.1 4.0 3.9 3.8 3.7 3.6 3.5 

Figure 2: Comparative ""H NMR spectra of pure (3-CD and a 1:1 
(3-CD-xylazine mixture showing (3-CD regions. 
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xylazine but the magnitude of downfield shifts was relatively 
small. These observations confirmed the formation of the P-CD- 
xylazine complex. 

Stoichiometry 

Of the various methods used to determine the stoichiometry and 
binding constant of host-guest complexes, NMR titration 
experiments are most common. These titration experiments give 
A5 data for many independent signals which can independently 
be used to determine the stoichiometry and binding constant. 
Several NMR versions of Benesi-Hildebrand equations [17] are 
used for this purpose and we have used Scott's method [18]. 
The Scott's equation gives the relationship between the 
apparent binding constant (K^), the complexation-induced shift 
at saturation (A5s), the guest or host concentrations (the species 
whose concentration is varied), and the observed shift change 
(A5obs)- Scott's equation for 1:1 P-CD-xylazine complex can be 
written as: 

[XZ] _ [XZ] 1 

NMR spectra of pure P-CD, pure xylazine and their mixtures 
having molar ratios (P-CD: xylazine) 1:0.25, 1:0.50, 1:0.75, 
1:1.00, 1:1.50 and 1:2.00 were recorded. NMR samples of 
mixtures were prepared by taking 9 mg of P-CD and adding a 
calculated amount of xylazine in 0.5 ml D2O. 

The A5h-3' aiid data was plotted in the form of [XZ]/ 

A5obs versus [XZ] (Figure 3) which gave linear fits confirming 
the 1:1 stoichiometry of the complex. The binding constant of 
the complex was calculated to be 86.45 M~^ 
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Figure 3: Scott's plot of the chemical shift changes of (B-CD cavity 


protons during titration with xylazine in D2O. 



^H-''H ROESY studies 

ROESY is an important tool for the study of large molecules 
and provides information on the through-space proximity of 
protons [19]. The host-guest interactions are displayed as inter- 
molecular peaks between cavity protons and part of the guest 
involved in complexation. The strength of the crosspeak is 
proportional to the inverse sixth power of the distance between 
the interacting nuclei, / oc 1/r^. However, for quantitative 
ROESY data analysis, one must understand the consequences of 
spin diffusion which occurs primarily for large molecules and 
long mixing times outside the "linear approximation" resulting 
in TOCSY artifacts. With a mixing time of around 0.2 s, 
TOCSY artifacts are generally not observed for large mole- 
cules (MW > 1200) but experiment requires several hours. On 
the other hand, ROESY spectrum of large molecules can be 
recorded in few minutes with long mixing time, e.g. 0.5s, but 
containing TOCSY artifacts. 

We recorded a ROESY spectrum of a 1:1 mixture of P-CD and 
xylazine with 0.5 s mixing time (Figure 4). The experiment was 
stopped as soon as the crosspeaks between xylazine and the 
cavity protons appeared. The region showing intermolecular 
peaks contained no TOCSY artifacts on one side of the diag- 
onal but artifacts interfered with intermolecular peaks on the 
other side of the diagonal. The intermolecular crosspeaks 
between H-1 and H-2 of xylazine and the cavity protons 
confirmed the presence of an aromatic ring inside the P-CD 
cavity. It was found that the heterocyclic ring was not involved 
in the complexation. Also, no interactions were observed 
between the methyl groups of xylazine and the cavity protons 
suggesting partial penetration of the aromatic ring but whether 
the aromatic ring approached the cavity from the wide or 
narrow side was not clear. Thus, two geometries for the P-CD- 
xylazine complex can be assumed (Figure 5). We established 
the structure of the complex using this ROESY data and molec- 
ular dynamics simulation studies in vacuum. Earlier, we estab- 
lished the structure of a fexofenadine-a-CD complex [20] on the 
basis of molecular mechanics and ROESY data recorded with 
0.5 s mixing time, but intermolecular peaks on both sides of the 
diagonal were clear in that case. 

Molecular dynamics 

Molecular dynamics simulations for only aromatic ring were 
performed. All the calculations were performed using CS 
ChemSD Pro (Cambridge Soft Corp.) in vacuum at 298 K. 
Initial coordinates for P-CD were obtained from the Cambridge 
databank. The published X-ray coordinates for hydrated P-CD 
[21] were used as starting point after removal of the water 
molecule coordinates. The structure of xylazine was minimized 
to a RMS value of 0.1 kcal mol~^A~^ using Allinger's force 
field. Simulations were performed by placing the aromatic ring 
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8.0 7.0 6.0 5.0 4.0 3.0 

Figure 4: Full ROESY spectrum of a (3-CD-xylazine 1:1 mixture showing intermolecular crosspeaks. 




Figure 5: Two probable modes of the inclusion of xylazine into the 
(3-CD cavity. 
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Figure 6: Coordinate system used to define the complexation process. 



of xylazine near the mouth of the P-CD cavity along the z-axis 
either on the narrow (N) or wide (W) side (Figure 6). 

every 10 iterations with 4000 steps of equilibration. Molecular 
The carbon skeleton of the P-CD was kept static but other atoms dynamics simulation results provided evidence that complexa- 
and the xylazine molecule were allowed to move. An iteration tion of P-CD and xylazine is energetically favored. In both the 
step of 1 fs was used and conformations were recorded after cases, the total potential energy of the complex was less than 
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the sum of the potential energies of the two components 
(Figure 7). 

The aromatic ring entered the cavity quickly and remained 
inside occupying various orientations by moving sideways, 
rotating or tilting, when xylazine was kept near the wider 
opening. The methyl groups always remained outside the 
cavity. In case of narrow side simulation, the aromatic ring also 
partially entered the cavity, with the methyl groups staying 
outside, but left the cavity after 2000 iterations and never 
returned. Figure 8 and Figure 9 show series of pivotal snap- 
shots of the two simulation trajectories. ROESY peak intensity 
ratios were then calculated from the interproton distances of 
two frames from narrow side (frame 105 ,1050 fs, 97.83 kcal/ 
mol; frame 118 , 1180 fs, 97.37 kcal/mol), and one frame from 
wide side (frame 132 , 1320 fs, 100.73 kcal/mol). 

Quantitative ''H-^'H ROESY analysis 

The ROESY peak intensities depend on the internuclear 
distances but are affected by several other factors and this is 
why the quantitative use of ROESY is generally avoided. Still 



there are numerous examples [22-24] where highly accurate 
structural information has been deduced by quantitative analysis 
of ROESY data. Macura and coworkers [25] and others [26] 
have shown that employing relative rather than absolute NOE 
intensities from within a given experiment can be used to calcu- 
late internuclear distances and vice versa with high accuracy 
using following equation. 



(2) 



^2 



where I\ and I2 are intensities of two ROESY peaks and r\ and 
r2 are distances between interacting protons in a given ROESY 
spectrum. 

To see whether this relation can be useful for the study of CD 
complexes, we performed a molecular dynamics simulation of 
the aspartame-P-CD complexation from the wider side. From 
the interproton distances, obtained from the lowest energy 
frame, we calculated the intensity ratios of all the intermolec- 
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2250 fs 3000 ts 4000 fs 

Figure 9: Snapshots showing the inclusion of xylazine into the p-CD cavity as obtained in the MD trajectory from narrow side. 



Beilstein J. Org. Chem. 2013, 9, 1917-1924. 



ular peaks. These peak intensity ratios matched quite well with 
those calculated from the reported ROESY spectrum of the 
mixture of aspartame and P-CD [27] but slight refinement of the 
lowest energy frame conformation gave peak intensity ratios 
which were in very good agreement. It was observed that peak 
intensity ratios for each interaction individually, for example, 
%-ortho-H-3'/%-ortho-H-5S OY summed intensity ratios (S/h-3'/ 
E/h-5') matched very well with experimental ROESY inten- 
sities. 

We then calculated the summed peak intensity ratio E/h-i/2^/h-2 
for two frames from the narrow side (Frame 105 and 118) and 
one frame from the wider side (Frame 132) obtained in molec- 
ular dynamics simulations. The interproton distances between 
aromatic protons of xylazine and cavity protons were obtained 
for the frames to be studied. All the H-H distances of each 
aromatic proton with seven H-3' and seven H-5' protons were 
obtained and their referenced peak intensity ratios (/h-i-h-3S 
%-l-H-5', %-2-H-3', %-2-H-5') were calculated using Equation 2 
which were then summed to give referenced S/h-i-h-3S 
I:/h-i-h-5', ^^h-2-H-3', ^^h-2-H-5'- The summation of I:/h-i-h-3' 
and S/h-i-h-5', as well as S/h-2-h-3', ^%-2-H-5' gave the total 
referenced /h-i and /h-2, respectively, from which the ratio was 
calculated and the results are given in Table 1 . The peak inten- 
sity ratio of H-1 and H-2 with cavity protons is closest for 
lowest energy frame from wider side suggesting that this must 
be the averaged ensemble conformation of the P-CD-xylazine 
complex (Figure 10). It must be mentioned that the potential 



energy for narrow side entry is lower than for the wider side and 
so it must be favored. But, unlike molecular dynamics simula- 
tions explicitly in water where contacts outside the cavity are 
self-compensative, the outside contacts in vacuum also 
contribute to the total energy and thus energy does not necessar- 
ily reflect the complexation energy. 



Table 1: Peak intensity ratios calculated from ROESY spectrum 
(experimental) and from conformations obtained by molecular 
dynamics simulations. 



Potential energy 
(kcal/mol) 



^H-1 / Ih-2 



Experimental 
Frame 105 (NS) 
Frame 118 (NS) 
Frame 132 (WS) 



97.83 
97.37 
100.73 



4.7 
8.7 
10.2 
5.0 



Conclusion 

The results of the study demonstrate that by a combination of 
quantitative ROESY analysis and molecular dynamics in 
vacuum, highly accurate structural information regarding the 
whole conformation of the inclusion complex can be obtained 
which otherwise only confirms the inclusion of the guest in the 
cavity. Studies on the determination of the absolute configur- 
ation of cyclodextrin complexes by quantitative use of ROESY 
data are in progress. 




Figure 10: Side and front views of the proposed conformation of the (3-CD-xylazine complex. 
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Supporting Information 

Supporting Information File 1 

The zip-archive contains an xls-file with the molecular 
dynamics wide side-energy-time plot, a molecular 
dynamics wide side trajectory in c3d format and the 
proposed beta-cyclodextrin-xylazine complex in c3d 
format. 

[http ://www.beilstein-j oumals . org/bj oc/content/ 
supplementary/1 860-5397-9-226-S 1 .zip] 
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